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The high statistics COMPASS results on diffractive dissociation ttN —► irnnN suggest that the 
isospin 1 = 1 spin-parity J PC = 1 ++ ai(1260) resonance could be split into two states: ai(1260) 
decaying into an S-wave pn system, and aj(1420) decaying into a P-wave /o(980)-7t system. We 
analyse the reaction by incorporating our previous treatment of resonant re-scattering corrections 
in the Drell-Deck forward production process. Our results show that the COMPASS results are fully 
consistent with the existence of a single axial-vector ai resonance. The characteristic structure of 
the production process, which differs in the two orbital angular momentum states, plays a crucial 
role in this determination. Provided the theoretical analysis of the reaction is done in a consistent 
manner, this single resonance produces two peaks at different locations in the two channels, with a 
rapid increase of the phase difference between their amplitudes arising mainly from the structure of 
the production process itself, and not from a dynamical resonance effect. In addition, this analysis 
clarifies questions related to the mass, width, and decay rates of the a\ resonance. 

PACS numbers: 12.40.Yx, 13.25.Jx 

I. INTRODUCTION 0”, 1 _ , and 0 + (the latter being somewhat ill- 

determined yet) form the ground-state of the light 
quark-antiquark spectrum. A newcomer in the 
family would be difficult to accommodate. 

2. It is peculiar to have two J PC = 1 ++ three-pion 
states, with identical quantum numbers, close in 
mass (within a full width of each other), with 
orthogonal decay modes, without the presence of 
some new quantum number. The Kg — K ° system 
led to decisive discoveries in fundamental physics; 
neutrino mixing is a spectacular current example. 
However, in the ai case, we see no candidate for a 
distinguishing quantum number. 

3. In the succesful quark-antiquark potential model 
approaches to the hadron spectrum, there is no sign 
that two bound states with identical quantum num¬ 
bers and comparable masses could be constructed 
from the u and d quarks (unless the potential is 
pathological). One would have to resort perhaps to 
a four-quark or molecule-type object. 

4. Regarding the COMPASS data, we may mention 
that a similar situation may exist in the case of 
the production of the 7 t 2 (1670) and a tt' 2 as s- and 
d-waves of 7 t/ o ( 1270). We defer this consideration 
and report here only on the a\. 

Our basic approach to high energy forward production 
of three pion states in pion-nucleon interactions is the 
Drell-Deck model 0. This model has been studied 
extensively in the production of the J p = 1 + pn system 
00, and here we extend the analysis to the J p = 1 + 
/o( 980 ) 7 t system. An important difference is that 
whereas the p7r system is in an orbital S-wave state, the 


Results of unprecedented precision are emerging from 
the COMPASS experiment at the CERN SPS where 
beams of 190 GeV pions interact with proton or nu¬ 
clear targets, providing new insight into the properties 
of light mesons 0. This investigation is of significant 
interest since it is bound to clarify low-energy hadron 
spectroscopy, where the situation of some states is some¬ 
what confused compared to expectations in a naive view 
of the standard model. 

In this paper, we concentrate on a question of impor¬ 
tance for many reasons: the isospin 1 axial-vector res¬ 
onance ai (reported in the Particle Data Compilation 
as a(l)(1260) [2j.) Evidence emerges in the COMPASS 
data for a new narrow J PC = 1 ++ axial-vector state 
with isospin 1, strongly coupled to the 7t/ 0 (980) system. 
This observation of a peak in the axial-vector two body 
7t/o(980) P-wave intensity at a mass of 1.42 GeV, com¬ 
bined with a phase motion close to 180° with respect to 
other waves, appears at face-value to mean that a sec¬ 
ond axial-vector resonance is present, close in mass to 
the known broad ai(1260) that couples mainly to the np 
meson channel 0. While these three features, i.e. two 
peaks at different masses and a rapid phase variation, are 
clearly present, there are reasons to be surprised, among 
which we mention the following: 

1. The ai(1260) is a central member of the axial- 
vector nonet, which, together with the J p = 
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/o(980)7r is in an orbital P-wave state, not studied to 
date. In addition, since we want to study resonances 
in these quasi two-body systems, we must modify the 
Deck mechanism with the proper corrections due to 
the re-scattering of these states. This is an inescapable 
physical consistency condition of the entire analysis. 

The unitary coupled channel approach that we devel¬ 
oped in the late 1970’s in Refs, [g], Q and [f|, should be 
an ideal way to show whether one resonance can produce 
mass peaks at different locations in the two decay chan¬ 
nels, along with a relative phase variation between the 
two channels, or whether the COMPASS data do require 
two nearby resonances with the same axial-vector quan¬ 
tum numbers in the three-pion system. In this paper, we 
demonstrate that the main features of the COMPASS 
data are compatible with a single resonance. We concen¬ 
trate on developing the theoretical method of analysis 
and its main consequences, leaving a detailed fit to the 
data to future discussion with the COMPASS collabora¬ 
tion. 

As we show, the double peak structure, by itself, 
requires a revision in the determination of the nominal 
second-sheet pole parameters of the a±, i.e. its mass 
and width. Owing to its large width, the a\ peak, 
as observed in various final states, appears distorted 
by several effects. In diffraction dissociation, as we 
have shown in Q, the structure of the Deck amplitude 
alters the resonance peak considerably. In r-decay 0, 
a distortion arises from phase space factors. Therefore, 
the apparent “mass”, as identified by the peak position, 
can vary considerably according to the production 
mechanism. As a byproduct of our present study, we are 
able to determine a new estimate of the branching ratio 
of the a\ into 7r/ o (980). 

In Sec. El we recall the basic facts about the Deck 
mechanism. We exhibit its P-wave structure which has 
not received much attention before now. We elicit a 
simple property that turns out to be the cause of the 
rapid rise up of the fan phase relative to pn. 

Since we are interested in the a\ which corresponds 
to resonant behavior of the pit and /o 7 r amplitudes, we 
must take into account the final state interactions be¬ 
tween these particles and/or the re-scattering corrections 
to the bare Deck amplitudes. It would be inconsistent 
physically not to perform such an analysis that incorpo¬ 
rates unitarity requirements. We published such work in 
the late 1970’s on the ai and on the AT1(1270)-A"1(1400) 
system, but our numerical results then were compared 
with data having much smaller statistics. In Sec. IIII1 
we describe the unitarization procedure, and we intro¬ 
duce the relevant physical parameters. We explain how 
we deal with a very restricted set of parameters when we 
introduce a single ai resonance. 

In Sec. HVI we present our results together with the val¬ 
ues of physical parameters involved, and we summarize 


our conclusions in Sec. E For the sake of clarity we col¬ 
lect the simple but necessary formulae in the Appendix. 


II. TWO-CHANNEL DECK AMPLITUDES 

We follow closely Refs. 0 , 0 and @. We consider 
the Deck amplitudes Td( 7 tN -a 7 T 7 T 7 rN) for production of 
the 7 T 7 T 7 r system at small momentum transfer and high 
incident energy (known as “diffractive production”). For 
the quasi-two body systems np and 7r/o, we denote 

Tp = T d (ttN -A 7 rpN) and 

= Td(7tN -A 7rf 0 N) . 

The reactions are represented in FigsJT](a) and (b). 




(b) 


FIG. 1: Deck production processes for (a) pn and (b) fon. 


The 7 xp case has been studied at length. Its amplitude 
is given by Eq. (2.1) of Ref. 0 as 

T& = gpnnKpfo) 2 l is 13 e btl a np ( 1 ) 

m*—t2 

where g p7r7r is the pnn coupling constant ( g 2 / ( 47 r) = 2.4), 
K p is the magnitude of the incident pion momentum 
in the p rest frame [sj, b is the slope of the ttN elastic 
diffraction peak, and a np is the np total cross-section. 
The invariants S 13 , <q and £2 are labeled in Fig.[T|(a). 

Similarly, the nfo production amplitude is 

If) — 9fomr 2 ~l~ ^136 lf ^7rp ■ (2) 

^7 r ^2 

This expression must be projected onto the orbital an¬ 
gular momentum l = 1 P-wave of the tt/q system for 
the purposes of our present investigation. Note here 
that the /o 7 T 7 t coupling constant, <?/ o7r7r , has the dimen¬ 
sion of a momentum. Choosing the average value of the 
/o —> 7 T 7 T width of 60 MeV, we obtain a numerical value 
<?/ 07r7r — 1-45 GeV. The other factors in Eq. ©, relative 
to Fig. [T] (b) have the same meaning as in Eq. |1]). 
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A. Background pn Deck amplitude 


fer. A convenient dimensionless expansion parameter is 


The Deck background amplitude has been well studied, 
where, by background, we mean the amplitude before any 
unitarization or rescattering correction is performed. We 
refer to our previous work in Ref. [gj and extract what 
is useful in the present analysis. We work in the final pir 
(or, more generally, resonance- 7 r) center of mass frame; 
M is the invariant mass of this system, and we denote 
by 9 the angle between the incident and outgoing pion 
momenta. 

In the limit of forward production (ti -A 0) and large 
s, the amplitude in Eq. © takes on a very simple form; 
see, e. g., Ref 0 

1 S 

T£ > (M 2 ,s,t 1 = 0)^a— - 7 , (3) 

M z — m z 

which is nothing but the (pir) S-wave projection Tj P _ 1+ 
of Eq.©, where a is a known constant that we discuss 
below. This S-wave projection has been used in previous 
calculations (e.g., Refs @] and 0 )- 

In the limit t\ —> 0 and large s, the Deck amplitudes 
for higher partial waves vanish identically. However, 
the J p = 1+ /o 7 r system is in an orbital P-wave. To 
address /o 7 r, we must extend the partial wave extraction 
calculations to finite values of t\ and s. The complete 
calculation of these amplitudes is presented in the 
Appendix. The important feature is that the higher 
partial wave amplitudes are of order t\/M 2 or M 2 / s 
compared to Eq. ©, and this S-wave amplitude is 
modified slightly. An immediate consequence is that 
fair P-wave production should have a noticeably smaller 
rate than the pn S-wave process, as is borne out in the 
complete calculation in the Appendix, and exhibited 
by the COMPASS data, where the intensity of the fan 
peak at 1.42 GeV, is lower than that of the pn peak at 
1.26 GeV by a factor of the order of a few ICC 3 . 

In the COMPASS experiment Ref 0 , the value of the 
square of the invariant total energy is s = 380 GeV 2 
while the momentum-transfer t\ in the smallest bin is 
t\ £ [—0.1,—0.13] GeV 2 (we follow the experimenters’ 
definition t\ = t — t m i n , except that, for our convenience, 
we work with negative momentum transfers). In this 
analysis, we are interested in values of M ~ 1 to 2 GeV. 
We notice that \ti\/M 2 M 2 /s and therefore the only 
relevant kinematic corrections come from the momentum 
transfer dependence (confirmed in the quantitative anal¬ 
ysis in the Appendix). We choose to work at the fixed 
value of t\ = —0.1 GeV 2 , and we checked that within the 
first t-bin (ii £ [—0.1,—0.13] GeV 2 ), our results do not 
vary appreciably. 


Bi = 


(M 2 — m 2 )' 


( 4 ) 


The J p = 1+ S-wave pit background amplitude is, to 
first order in 0 1; 


T Deck 
S 


(M 2 — m 2 ) 


2 ) e p,A^ 1 + u 


_ 1 (3 M 2 +m n , _ 

2 UU (M 2 -m 2 ) 1 -y 


( 5 ) 


where E n and E p are the pion and p energies in the p7r 
rest frame, and where 


y=p n /E n , 


( 6 ) 


is the pir phase space factor, being the pion momen¬ 
tum in the pir rest frame. 

The J p = 1+ P-wave /o7r amplitude is, at the same 
order in 0i 


T, 


Deck 


1, 


3 s 

01 X 

2 (M 2 - ml) 

( (3 M 2 + ml) 

^ I 

\ (M 2 - ml) 

E-k 


,—2 1 

(-k 2 

y y 2 


ln(^))(7) 

i -y 


where E^, Ef 0 are the pion and /o energies, p n the 
pion momentum in the fair rest frame and, as above, 

y = Pir/En. 



FIG. 2: Behavior of the background P-wave fon Deck ampli¬ 
tude above threshold, showing the zero of this amplitude near 
1.38 GeV. 


B. pn and /o7r background amplitudes to first order 

In the Appendix we derive expansions of the back¬ 
ground amplitudes to first order in the momentum trans- 


This latter expression contains a major clue to our in¬ 
vestigation. Of course, we shall unitarize these formulas, 
and combine them properly, but the simple expression of 
this /q7t Deck amplitude explains an important feature 
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of the data. Indeed, the right hand side contains the 
factor (3 M 2 + ml)/(M 2 — ml) — Ef 0 /E v , plotted in 
Fig. 0. This factor is negative at low values of M 
(since mf 0 /m n > 3), but it vanishes at some point 
near M ~ 1.38 GeV and becomes positive afterward. 
Furthermore, if we give this term some small imaginary 
part, its phase will suddenly switch from —180° to 
zero. This sudden and rapid phase variation is not a 
dynamical effect in the sense of a resonant phase, but 
it originates in the structure of the dynamical process 
by which the / 0 7 r state is produced. A similar term 
is present also in the S-wave, Eq. ©, but there it is 
completely hidden by the dominant term. 

Another interesting qualitative feature of Eq. 0 is 
that it grows in the region of interest (M ~ 1.2 to 
1.4 GeV) and therefore tends to push a resonance peak 
upward in M, and, because of the zero, to produce a 
(non-resonant) peak or maximum between threshold and 
M = 1.38 GeV. 


C. Parameters and normalization 


Keeping in mind the parameters introduced in Eqs. 0 
and ©, our two J PC = 1 ++ amplitudes are 

f T Deck (pn) \ _ 2iy/2sN / 9pirnKp(J 7r pTp 7r \ 

\ T Deck(fo^) ) ( M 2 - ml) V 5/ o7r7r CTirp Tfo-K ) ' 

( 8 ) 

where, T pn and T/ o7r can be read off from Eqs. © and 
©. The structure remains the same after we unitarize. 
The normalization factor N is irrelevant and is taken 
equal to 1 here. In a more complete analysis N will 
ensure precise quantitative fits to the data. 


After integration over phase space, the normalized dif¬ 
ferential cross sections, expressed in GeV units, are 


d<7 , _ 1 1 Qa* 7r 

dtidM ltl ~ 2 0.3893 (4t r) 3 s 2 


\T D eck{a*n)\ 2 


(9) 


where a* stands for p or /o, and (] a *n is the magnitude 
of the a* in the a*n rest frame. 


III. UNITARIZATION 

For all theoretical and technical details about multi¬ 
channel final state unitarization, we refer to the litera¬ 
ture, in particular to Ref. 0 where the general analysis 
is done thoroughly and to Ref. @ where a specific ap¬ 
plication is made. Let us simply recall that if S is the 
(two-channel) strong interaction S matrix, then the mili¬ 
tarized Deck amplitude 7z>, which we can write as a two 
dimensional vector as in Eq.©, has a right hand unitar- 
ity cut along which it satisfies the relation 

Tv + = STv~ ( 10 ) 


7r> + and 7x> being the values of the unitarized Deck 
amplitude above and below the cut. 


A. The resonant pn and fon systems 


In order to unitarize the Deck amplitude by taking 
into account resonant inelastic (or coupled channel) final 
state interactions between the pn and fpn states, we in¬ 
troduce Chew-Mandelstam functions [ll[ C p and C ?. By 
definition these are analytic functions of the invariant en¬ 
ergy squared s of two particles, with a right hand branch 
cut where the imaginary part is equal to the phase space 
factor 2p/y/s, p being the c.m. momentum: 


G mj ^(s) — C(s, m, p) — 

— — { - -[(to + p) 2 - s] 1 / 2 [(m - p) 2 - s] 1 / 2 

7T S 

[(to + p) 2 — s] 1 / 2 + [(to — p) 2 — s] 1 / 2 
2 (mp) 1 / 2 

to 2 — p 2 to to 2 + mu 2 to 1 

2s p 2 (to 2 — mu 2 ) p 2 J 


( 11 ) 


The function C P (M 2 ) (where M 2 is the square of the 
c.m. energy) is written at length in Eq. (3.8) of Ref. 

is given by a slightly different formula since nfo is in 
a P-wave. 


The P-wave nature has a few simple consequences. An 
elastic nfo —> nfo amplitude behaves, at threshold, as 

_2 rn ^ W 2 ~ ( m f +m^) 2 ] [. M 2 - (to/ - to w ) 2 ] 

- > 

( 12 ) 

and an inelastic amplitude such as np —> nfo behaves as 
p n f 0 (M 2 ). This last property can be seen in the Deck 
production amplitude Eq. © which is proportional to 
y = 2p/M as y —> 0. 

In practice, since Eq. ED vanishes at M 2 = 0 , and 
therefore p 2 / o (M 2 ) x C(M 2 \ mf, m^) is regular and an¬ 
alytic in the cut plane from (to/ + to w ) 2 to infinity with 
a discontinuity equal to 2p 3 /M , it is more convenient for 
us to work with the two Chew-Mandelstam functions of 
different dimensions: 

Ci(M 2 ) = C P (M 2 ) = C(M 2 ; to p , TO-n-) , (13) 


and 

C 2 (M 2 ) = C f (M 2 )=pl fo (M 2 )C(M 2 -m f ,m n ) . 

(14) 

These have different dimensions, but they will appear 
more symmetrically in the formulas, making these easier 
to read. They are also associated with constants of dif¬ 
ferent dimensions, and we take care, inside the complete 
formulas, to check the presence of the appropriate pow¬ 
ers of the momenta. This is an alternative to the more 
general treatment described at length in the literature, 
for instance in Eqs. (3.1) to (3.4) of Ref. 0 . 









5 


B. Unitarization procedure 


We follow the extensive coupled channel unitariza¬ 
tion procedure presented in Sec. Ill of Ref. [gj. We 
assume that the p and the fo are approximately stable. 
Therefore C p and C? are defined analytically as above. 
Furthermore, we assume that the fo is an elastic 7T7 t 
resonance. We know that this is not true. There is a 
noticeable KK inelasticity, but it is not predominant. If 
it turns out to be necessary to incorporate inelasticity, 
the simplest procedure consists in introducing the 
KK partial decay mode through a third channel I\K , 
with the corresponding Chew-Mandelstam function C K } 


Our basic assumption, in the present analysis, is 
that there is a single ai resonance whose (unique) 
second-sheet pole parameters we shall determine. Since 
we are dealing with a two-channel case, we parametrize 
the coupled pit and foit final state interactions (or 
rescattering) via this resonance. 

In order to do this, we introduce a K matrix, as in 
Eq.(3.14) of @ 

( 9l 9192 \ 

XT ) • (15) 

si — AI ‘ 1 2 s-L-M ' 2 ) 

As mentioned above, since the two channels have the 
same total angular momentum but different orbital 
angular momenta, the a± —> pit coupling constant g\ has 
the dimension of an energy, while the a\ —>• foit coupling 
constant g 2 is dimensionless. The squared mass si will 
be related to the a\ resonance mass. 2 All our units 
are in GeV. We work with the parameter 7 = < 72/51 hi 
discussing the results. 


The crucial tool to treat coupled channel final state 
interactions is the ZTmatrix, which, in this case, is given 
explicitly in Eq.(3.15) of @ by 


D{M 2 ) 


1 

V 0 (A/ 2 ) 


gi -52(^1 - M 2 - a 2 C 2 ) A 
52 51 (si - M 2 - a 2 Ci J 


(16) 


1 We actually built an unstable Chew-Mandelstam function for the 
fo following what we did in Ref. [?J in the case of the K 1(1270) 
and K1 (1400) strange axial vector mesons. This procedure re¬ 
produces the decay width of the fo as well as the observed n — n 
inelasticity El at the opening of the KK channel, in analogy 
with the empirical Flatte formula El- However, we found very 
small deviations in the results compared to the stable case be¬ 
cause the coupling of the ai resonance to the /o7r channel is 
small. We might restore inelasticity in a subsequent analysis. 

2 One can perfectly well operate in a more standard way, but the 
momentum factors generate slightly complicated formulas. Here, 
one can include the p factors due to the n — fo amplitude directly 
into the constant g 2 which would become an expression of the 
type pg 2 - This would ensure the correct P-wave behavior of the 
7r — fo channel, without having to change the Cf function. 


where a 2 = gf + gf > an d where the energy denominator 
function V 0 {M 2 ) is 

Vo (M 2 ) = ( Sl - M 2 - glC^M 2 ) - g 2 2 C 2 {M 2 )) . (17) 


The function Vq (M 2 ) is of fundamental importance. 
It contains all the information (that we put in) on the 
coupled-channel pit-foit strong interaction. Equation 
(14) of Ref. 0| shows how the matrix D gives direct 
access to the strong interaction t— matrix (whose first 
diagonal element is e lS sin S in the elastic region). 
Vo (M 2 ) is an analytic function which possesses the pit 
and foit branch cuts from [m p + m^} 2 to infinity and 
from [ mf 0 + m,r] 2 to infinity. Its second sheet pole 
determines the nominal position and width of the a\ 
resonance. The function 1 ? 0 (M 2 ) is illustrated in Fig. ([3]). 



FIG. 3: Behavior of 1 /|77>o (A/ 2 )| 2 as a function of the energy 
M, in the case of an input second sheet pole corresponding to 
m ai = 1.4 GeV and r ai = 0.3 GeV as used in this analysis. 


As in 


we define the unitarized Deck amplitude as 


Tfi){M 2 ) = T d (M 2 ) - —D(M 2 ) x 


ds 


(m p +m ^) 2 


, ImD{s')T D (s') 
{s' - M 2 ) 


• ( 18 ) 


Here T^{M 2 ) is a two-dimensional vector and Tjj{M 2 ) 
is the “background” Deck amplitude discussed above. 

Equation (1181) is the Deck amplitude with rescattering 
corrections taken into account. It has the same left-hand 
singularities as Td(M 2 ) (for instance the pole at M 2 = 
m 2 ), it satisfies the unitarity relation Eq. (flOl) . It reduces 
to Tjj{M 2 ) if the strong interaction amplitudes vanish. 
In the Appendix of Ref. Q we justify this interpretation 
qualitatively within a simple effective Lagrangian model. 
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C. Direct production amplitude 

There is one more ingredient. In the same effective La- 
grangian model mentioned above, we show that in most 
cases, such as the present problem, one must also incor¬ 
porate a direct production contribution of the form 

T(M 2 ) = T£(M 2 ) ± D(M 2 )P(M 2 ) , (19) 

where P(M 2 ) is a constant or first degree polynomial. 
The corresponding process is illustrated in Fig. (0. 



FIG. 4: Diagram representing direct production of ai reso¬ 
nance through coupling to the diffractive exchange mechanism 
(pomeron. 

There are various arguments in favor of this term, such 
as high energy behavior. The most intuitive is based 
on an examination of a one channel case where the ex¬ 
plicit result, the Omnes-Mushkhelishvili equation 0, 
has been known for a long time. If the resonance is 
not a dynamical effect of the forces, but rather a CDD 
pole, i.e. a system with appropriate quantum numbers 
whose origin comes from another system, such as a quark- 
antiquark state, then a form of the type of Eq. m will 
produce a unitarized amplitude behaving as e lS cos <5 with 
a zero near where we expect the resonance to appear 
(<5 ~ 90°), whereas the direct production term behaves 
as e I<5 sin 5 with a bona fide peak. A combination of the 
two terms is necessary in general. When we display the 
various stages of before and after unitarization, we see 
how this zero appears, why direct production is necessary 
by itself, and also why both contributions are necessary 
in order to build the observed resonance peak. 

In the present analysis, the direct production term 
must also be diffractive and obey the unitarity con¬ 
straints. We choose an amplitude of the form 

T dir (s,M 2 ) = ^ay(£), (20) 

where G is an arbitrary constant. Ideally, one would ex¬ 
pect the two constants which appear here /i and ft to 
be proportional respectively to g\ and g 2 , which are the 
couplings of the ai to the pn and fon channels. This 
proportionality turns out not to be possible, although 
the values we find have appropriate orders of magnitude 
(i.e. I/2LC |/iI)- However, as we have pointed out pre¬ 
viously [7|, there is arbitrariness in the definition of the 
direct production amplitude, and it is not of primary im¬ 
portance for our main conclusions. 


IV. ANALYSIS AND RESULTS 

In tis section we present the main results of our 
investigation. We have not attempted a detailed com¬ 
parison with data, based, for example on a minimization 
procedure. Nevertheless, there are some salient points 
that can be made. 

Our first ambition has been to select appropriate values 
of the ai parameters that provide a good global descrip¬ 
tion. This goal was not as trivial as we thought. In our 
previous paper on the ai @ we had been able to explore 
a wide range of values of the ai mass and width. It turns 
out that the COMPASS results fix these parameters in a 
more stringent way than when we dealt only with the S- 
wave pn system, be it or not associated with other S-wave 
channels (such as K*K). Here, the acceptable mass and 
width of the ai, as defined by the position of the second 
sheet pole , which are the actual parameters of our calcu¬ 
lation (i. e., after fixing them, we determine the values 
of Si and <71 in Eq. tfTSlP . turn out to be quite restricted. 
Despite the fact that it is more oriented towards qualita¬ 
tive aspects of the COMPASS data than to quantitative 
fits, our analysis indicates that: 

M(ai) ~ 1.40 ±0.02 GeV 
r(ai) ~ 0.30 ±0.05 GeV. 

These values of mass and width correspond to values of 
the parameters Si ~ 2.002 GeV 2 and g\ ~ 0.732 GeV. 
We have varied the parameters in order to understand 
the global trends, but in this first calculation we prefer 
to keep them fixed: M(a\) = 1.4 GeV, T(ai) = 0.3 GeV 
(within the uncertainty ranges noted above). 

The interesting parameter to vary is the ratio 7 = 
52/51 in order to find the range of values that produce 
two peaks with appropriate characteristics: the fon peak 
occurs at higher mass than the pn peak, and the ratio of 
maximal intensities of these peaks, i.e. pn/ fon, falls be¬ 
tween 1,000 and 500, as indicated by the available data. 
These requirements lead to negative values of 7 in the 
range [—0.1, —0.055]. In other words, the ai couplings to 
pn and to /o7r have opposite signs! The fact that 7 is neg¬ 
ative is quite easy to understand. If we look at Fig. m, 
we see that the relative sign of the P-wave and S-wave 
amplitudes determines whether the phase switches from 
— 180° to 0 , i.e. upward, or from 0 to —180°, i.e. down¬ 
ward. For the same reason, the P-wave peak will tend to 
move to higher energies than the S-wave peak. In other 
words, the data impose that the relative couplings have 
opposite signs. We choose as our central value 

7 = 52/51 = -0.08 , 

and we show how the energy variation of the phase 
changes with this parameter. 

The next step is to determine the amount of direct 
production necessary to fix the two peaks, one in pn, 
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the other in fon, at their desired positions, i.e. M = 
1.26 GeV for pn and M = 1.42 GeV for f 0 n. Values such 
as 

Go np fi = 120 and Go 7Tp f 2 = 5.5 

in Eq. O ensure good positions for the two peaks, and 
this situation is stable when one varies the parameter 7. 
Two facts seem important. First, the parameters which 
are essential in fixing accurately the positions of the 
two peaks are the direct production amplitudes above. 
It is not possible to change appreciably the position of 
either of the peaks by varying the parameter 7. Some 
variation is possible if one changes the mass and width 
of the ai, but this approach leads to a deviation from 
the desired values of 1.26 and 1.42 respectively. The 
second point is that the ratio of direct production to 
the background Deck amplitude is consistent with the 
value we obtained some time ago [6] in our analysis of 
data at much lower energies, with smaller statistics (see 
Ref. [H[). In our opinion, this is a positive outcome. 
In fact, in the present analysis, the inclusion of the 
fon channel is a small perturbation on the dominant 
S-wave pn amplitude, and it is satisfactory to see that 
the results confirm this. 

It is interesting to see how the two peaks build up, 
now that we know in what range we can choose the 
parameters. We plot in Fig. © the shapes of the rho 
intensity as a function of the energy M, for the various 
stages of the calculation. We notice that the pure Deck 



FIG. 5: Evolution of the pn differential cross section as a 
function of the energy M, in three cases: Background Deck, 
Unitarized Deck as in Ea. dlSI) . and final version including 
direct production. 


background does not produce a resonant shape. The 
unitarized amplitude shows effectively the ~ cos S zero 
we mentioned above in the one-channel case (to which 
this problem is actually very close). Finally, the direct 
production produces the observed peak, at the observed 


position. 

A similar set of curves for the fan intensity is shown 
in Fig. ©. Notice that the form of the Deck background 
by itself appears to simulate a narrow resonance peak at 
threshold (of course without any accompanying phase). 



FIG. 6: Evolution of the fon differential cross section as a 
function of the energy M, in three cases: Background Deck, 
Unitarized Deck as in Ea. dlSI) . and final version including 
direct production. 

Because the fan peak has a small intensity, in order 
to see properly the separation of the S-wave and P-wave 
peaks, we prefer to multiply the fon intensity by an ap¬ 
propriate coefficient of 650 such that it appears on an 
equal footing with the pn peak. This result is shown in 
Fig. 0. 



FIG. 7: pn and fon differential cross sections as a function 
of mass. The second is multiplied by a factor of 650 to make 
the figure easier to understand. 

The separation of the positions of the two peaks is 
evident. The width of ai(1260) is about twice the width 
of ai(1420) in the calculation, in agreement with the 
COMPASS observation. The ai(1420) peak is also more 
symmetrical, with width about 0.14 GeV. One inter¬ 
esting detail is that the lower end of the fon intensity 
exhibits a (tiny) peak at around 1.2 GeV owing to the 
zero emphasized in Fig. (0, which in turn is part of the 
dynamical structure of the P-wave Deck amplitude that 

































is responsible for the sharp rise of the phase difference 
between the ffi r and pit amplitudes. Perhaps this small 
peak near fait threshold is accessible in the experiment. 

We display in Fig. ([5]) a set of phase differences 
between the ffit and pit amplitudes for three values of 
our parameter 7 . The value 7 = —0.08 is, in some sense, 
our “central” value (subject to more refined analyses). 
We also plot the phase for a smaller value 7 = —0.95 and 



FIG. 8 : Three phase-differences between foit and pit, in the 
1.2 to 1.6 GeV region. The medium curve corresponds to 
our “central” solution 7 = —0.08 the two others (marked) to 
other values in the range of interest. In all three cases, the 
other results of the calculation remain practically unchanged: 
peak positions and intensities. 


for a larger value 7 = —0.065. For all three, we observe 
the spectacular rise of the phase due to the structure of 
the P-wave Deck amplitude. Notice that the larger value 
7 = —0.065 seems to be close to a limit beyond which 
our calculations develop unwanted effects. On the other 
side, for smaller values 7 =< — 0 . 1 , the phase tends to 
flatten and to be less dramatic in the transition region. 
Curiously, for all sensible values of 7 , all curves seem to 
cross at the same point near 65° at a mass of ~ 1.25 GeV. 

Empirically, it may be that the rise of the phase 
occurs somewhat too soon as compared to the data. On 
the other hand, ours is a first order calculation of the 
P-wave foit amplitude. The expansion parameter is of 
the order of 0.1 so that we can expect deviations of at 
least 10 % in further calculations, which we can perform 
as needed. 

Finally, an analysis such as this one can pin down the 
branching ratio for the ai decay into ffit. A precise value 
would come out of a more quantitative fit to data. It 
depends strongly on the factor 7 , and, for the time being, 
we prefer to state that the branching ratio (as compared 
to the dominant decay mode a\ —> pit) is of the order of 
10 " 3 . 


V. CONCLUDING REMARKS 

Our initial motivation for this study came from the 
fact that we found it implausible that two ai resonances 
could be so close in mass. This question led us to a 
new investigation of coupled channel rescattering effects 
in the Deck diffractive dissociation model. We were 
fortunate to rediscover in this example that, although 
a peak is often related to a resonance, its exact posi¬ 
tion depends also on the dynamics of the production 
mechanism by which it was created. The present case 
is interesting since it is the same diffraction dissociation 
process that is at work, but with radically different 
behaviors in different orbital angular momentum states 
of the dissociated products. Now that the calculation 
has been performed, this result seems quite reasonable, 
if not obvious, but we did not expect such important 
quantitative effects as a 160 MeV mass difference be¬ 
tween two peaks that correspond to the same resonance. 

One might question whether it would nevertheless be 
possible to have, within this theoretical framework, two 
independent J p = 1 ++ resonances so close to each other. 
At the onset of this project we did try the hypothesis 
of two resonances. But while the two-peak structure is 
easy to obtain, the phase difference is not. We had not 
yet understood fully the mechanism of the P-wave in 
the Deck model. So, in order to make this model cope 
with two resonances and agree with the COMPASS data 
one would have to be extremely lucky. First these two 
resonances should both be strictly coupled to different 
channels: one to pit and the other to ffit. In other 
words, some strict selection rule should be at work, 
otherwise the P-wave effect, that in the present case has 
been able to “double” a resonance peak might also act 
for each individual peak and predict not two peaks but 
three or four! 

As a last remark, we have analyzed the high energy 
production of the ai in a strong interaction process. If, 
idealistically, one were to perform a low energy itp and 
7 r/o elastic scattering experiment, one would observe a 
single resonance peak as a function of the center-of-mass 
energy, i.e. that of Fig. ©, at a position of, e.g., M=1.36 
GeV with width T = 0.31 GeV. 
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Appendix A: Kinematics and the partial wave 
decomposition of the Deck amplitude 

The process under consideration is ap —> a*nN. Four 
momenta of initial and final state particles are labeled p 
and q respectively. Three momenta are denoted by bold 
characters p. Norms of the three momenta are denoted 
by normal characters: p and q when no confusion can 
arise. 

We use A for the system (a*7r) and M for its invariant 
mass, which we also refer to as M A when it it useful. The 
invariants of interest are: 



a p —a* n N 


S (Pa H" Pp) (^a “1“ Q.TZ QN ) 

, (Al) 

SttN — {Qn + Qn) 2 = i 

(A2) 

Sira* = (Qit + Qa *) = -^A ? 

(A3) 

s a *N = (qa* + (In) 2 , 

(A4) 

taa* — (Qa* Pa) j 

(AS) 

t\ = tpjsj = (qN Pp) 5 

(A6) 

tan — (qn Pa) 

(A7) 

Note that 


Sna* t aa * “h t a 7T — Tfl> a "t" Tn a * "h ^n tpN 5 

(A8) 

and 


&na* + s nN + $a*N = S TTl + 171^ 

(A9) 


t CHANNEL 

FIG. 9: t-channel polar angles. 

Expressed in terms of polar angles, Eq. (A.5) becomes 

taa* =m 2 a ,+m 2 a -2E a *E a + 2\q a *\\p a \cos6 t . (A17) 

The relation of the angle \ between p a and qN, and the 
polar angles represented in Fig. ©, is 

cos ip = cos 0 t cos x + sin 9 t sin y cos cp t , (A18) 

where ip is the angle between q a » and qN- 


1. t-channel Variables; Rest FYame of A 


2. Amplitudes in the Deck model 


In this calculation, the interesting variables are the set: 
s, and t p N = t\ and the t-channel scattering angles, 
defined in Fig. [9l in the rest frame of A where q^ = 
q a * + Qtt = 0. In terms of these variables, the energies 
and magnitudes of the three-vector momenta of a, TV, a* 


and 7T are: 

E a = (M A + rr? a — ti)/2M A , (A10) 

Pa = \ l/2 {M 2 A ,m 2 a ,ti)/‘2M A = |p a | , (All) 

E n = (s — m 2 N — M A )/2M A , (A12) 

q N = X 1/2 (s,m 2 N ,M A )/2M A =\q N \ , (A13) 

E a - = (s - m 2 N - ml)/2M A , (A14) 

IQttI = |q a *| = A 1/2 (M^,to 2 ,,to 2 )/2 M a , (A15) 

E n = M A -E a , , (A16) 


where 

A (A 2 ,S 2 ,C 2 ) = [{A + B) 2 — C 2 ][{A — B) 2 — C 2 ] 


As can be inferred from Eq. JT]) the quantities of inter¬ 
est in this analysis originate from the ratio 


SnN 

h ~ ml 


(A19) 


where now we specialize to the case of an incident pion 
(to q = m,r))- Our aim is to obtain an expansion of this 
amplitude to first order in M 2 /s and in t\/M 2 . The 
following dimensionless quantity occurs repeatedly: 


Bi 


ti 

(M 2 — ml) 


(A20) 


The following expansions are useful: 


E 


a 


Pa 


(A/ 2 + ml- ti) 

2 M 

M 2 ~ ml f _ (M 2 Ami) \ 

2 M V (M 2 -m=) V' 


(A21) 
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e n 


<lN 


_s_ _ M 2 + m 2 N 

2 M [ s h 

s , _ M 2 + m 2 N 
2 s ’ 


(A22) 


Consider the numerator of Eq. IA19I We have: 


s-kN (.QttEqn') m -\-m2 7r -\-2E 1T E2q 7T qjj cosip , 

(A23) 

where ip is the angle between qN and q„-, as given by 
Eq. (1A1SI) . 

To calculate the angle x between qN and p a as defined 
by Eq. (fAISlt . one can make use of the identity: 

(QN ~ Pa ) 2 + {QN + 9o* + Qir) 2 + {QN ~ Pp ) 2 = 

(qN - Pa ) 2 + S + ti = m 2 N + m 2 + m 2 + M 2 (A24) 


One obtains, to first order in both ti/M 2 and M 2 /s, 


cos x = 



2 M 2 

(M 2 — m 2 ) 


(M 2 + m 2 N ) ^ 


(A25) 

where the quantity ©i has been defined in Eq. (IA20I) . 
Since we perform this calculation at the experimen¬ 
tal lower value of the momentum transfer, i.e. t\ = 
—0.1 GeV, whereas the experiment is performed at 
s = 380 GeV 2 , the last term on the right is negligible 
compared to the second one, and this feature remains 
true in all the calculation. We need to consider only 
the first order terms in ti/M 2 . It is understandable that 
cos x should be close to —1. since we consider the forward 
direction in high energy production. 

Using Eq. (1A18I) . we therefore obtain 


S-kN = JJ ( K - Pn COS 0(1 + 01 


2 M 2 


(M 2 — to 2 ) 


) (A26) 


Notice that we have not taken into account the second 
term on the right of Eq. (IA18I) since it does not contribute 
to the partial waves we are interested in (we consider 
only the m = 0 amplitudes). 


The denominator of Eq. (IA19I) is simpler. One obtains 
from Eq. (lAlTl) 

t 2 -m 2 = (p a -q a *) 2 -rn 2 = m 2 a *-2E a E a ,+2p a p 7I cosd 

(A27) 

The presence of a non vanishing value of t\ can be read 
off from Eq. (IA21I) . One obtains: 


t 2 - 


M 2 — m 2 


(K ~ ©1 E a , ) - p v COS 0(1 - 01 


(M 2 + TO 2 ) 
(M 2 — TO 2 ) 


|A28) 


Altogether, the Deck amplitude of interest Eq. (IA19I) 


to first order in t\ is 

rjiDeck __ S nN _ _ S _ 

t 2 — m 2 ( M 2 — m 2 ) 

/ E n — pn cos 0(1 + 0i )) V 

^ {E„ - 0i E a .) - p„ cos 0(1 - 0i ) f 

In the limit |ti| = 0 this expression reduces to 
— s/(M 2 — m 2 )) which is a pure S-wave. All other par¬ 
tial waves vanish in the limits s —> oo and t\ —> 0. Since 
Eq. (IA29I) is a rational function of cos0, the calculation 
of its partial wave projections is straightforward. 


3. Partial wave projections 


For convenience, we set: 


a 

b 

V 


l + 0i 


l-0i 


2 M 2 

{M 2 - ml) 
(M 2 + m 2 ) 

(M 2 — m 2 ) ’ 


Pnb 

E 1 r — 01 E a * 


and 


(A30) 


The S-wave projection of Eq. (IA29I) is obtained by in¬ 
tegration over cos0: 


rjiDeck _ _____ 

s “2 (M 2 - ml) 

( 2a 0 i ( F (3M 2 +m 2 ) a \ 

\b E„- 0i£ a . V ^ b(M 2 - ml) b a ) 


(— In 

y 



(A31) 


With the same notation, the P-wave projection of 
Eq. (IA29I) is obtained after multiplying by cos0 and in¬ 
tegrating over cos0: 

rjiDeck 0 _^_ 

p 2 (M 2 -ml) 

01 / (3M 2 + to 2 ) a \ 

(E v -e i£ fl -) 1 *b(M 2 -ml) b a *J 

(— + ^h\ ±1 )) ■ (A32) 

y y z 1 -y 

Notice that the above expressions are strictly the S- 
and P-wave projections of Eq. (IA29I) . The terms of order 
(0i) 2 are incomplete since terms of that order would also 
occur in Eq. (1A29I) . Therefore, to first order in 0! the 
S-wave amplitude is: 


rjiDeck __^_ 

S ~~(M 2 -m 2 ) 

/ 0i/ (3M 2 +m 2 ) E p \ 

\ 2 \ (M 2 — ml) E«) 


(Iln(l±l))W) 

y i -y / 
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and y = p v /E n . 


The P-wave amplitude is 

rjiDeck _ ^ 


2 (M 2 — to 2 ) 


( (3 M 2 + ml) E 


V (M 2 - ml) E, 




01 X 


-2 


y y‘ 


A + y 

i -y 


)).(A34) 
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